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Abstract 

We consider lepton flavour violating decays of the top quark, mediated by four-fermion operators. We compile 
constraints on a complete set of SU(3) xU(l)-invariant operators, arising from their loop contributions to rare decays 
and from HERA’s single top search. The bounds on e-/r flavour change are more restrictive than i-r; nonetheless 
the top could decay to a jet -l-e/2 with a branching ratio of order 10“^. We estimate that the currently available 
LHC data (20 fb~^ at 8 TeV) could be sensitive to BR{t —>■ efi+ jet) ^ 6 x 10“^, and extrapolate that 100 fb“^ at 
13 TeV could reach a sensitivity of ~ 1 x 10“®. 


1 Introduction 

Lepton Flavour Violation(LFV) [1], meaning local interactions that change the flavour of charged leptons, should 
occur because neutrinos have mass and mix. This motivates sensitive searches for processes such as ^ ey [2] 
and /i—e conversion [3]. However, the model responsible for neutrino masses is unknown, so it is interesting to 
parametrise LEV with contact interactions, and to look for it everywhere. In this context, the LHC could have the 
best sensitivity to LEV processes involving a heavy leg, such as the Z [4-6], the Higgs [6-9], or a top [10]. In this 
paper, we study the LEV top decays t —>■ qe^yP^, where q G {u,c}. 

We suppose that these decays are mediated by a four-fermion interaction, and outline in Section 2 the current 
bounds on LEV branching ratios of the top. The bounds arise from rare decays and HERA’s single top search, and 
are discussed in more detail in the Appendices. We find that, while these bounds place strong constraints on some 
specihc Lorentz structures for the 4-fermion interactions, they still allow for t —>■ qe^fjP decays with rates within 
the LHC reach. In Section 3, we estimate the LHC sensitivity to t —>■ qe^fjP, with 20 fb“^ of LHC data at 8 TeV. 
This estimate relies on simulations of the background and signal, and is inspired by the CMS search for t —>■ Zq [11]. 
The extrapolation to higher energies and luminosities is discussed in Section 4. 

Qitarfc-flavour-changing top decays, such as t —>■ cZ and t —>■ he, have been studied in the context of explicit 
models [9,12,13] or described by contact interaction parametrisations [14], and have been searched for at the 
LHC [11,15-17]. Quark-flavour-changing (but lepton-flavour-conserving) three-body decays of the top, t —^ off, 
where / is a lepton or quark, have also been calculated in explicit models [18]. Top interactions that change quark 
and lepton flavour, and, in addition, baryon and lepton number, have been explored in [19] and searched for by 
CMS [20]. In models with weak-scale neutrinos N [21], there can be lepton number- and flavour-changing W decays: 
W~ —>■ N£ —>• qq'i'i, which could appear in the final state of top decays. In the presence of this decay, ti production 
could give a final state with three leptons, missing energy and jets, as in the decay we study (see Figure 2). However, 
a different combination of leptons and jets should reconstruct to the top mass. Finally, Fernandez, Pagliarone, 
Ramirez-Zavaleta and Toscano [10] studied almost the same process as us, t qr^yp, but mediated by a (pseudo)- 
scalar boson. They obtain separately the low-energy bounds on the quark-and lepton-flavour-changing couplings of 
their boson, and obtain that LEV top branching ratios can be ~ 10“® if the boson mass is ^ 2mw- 
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2 Current bounds 


We are interested in the decays of a top (or anti-top) to a jet and a pair of oppositely charged leptons of different 
flavour. In this work, we focus on the processes t —>■ where q G {u, c}, because the e and /r are easy to 

identify at the LHC, and e O /r flavour violation is the most strictly constrained at low energy. We leave the decays 
to qe^T^ and qqt^r^ for a later analysis. 

We suppose that these decays are mediated by four fermion contact interactions. A complete list of the SU{3) x 
U{1) invariant operators that we study is given in Appendix A. We do not impose SU(2) on our operators, because 
the scale we will probe is not far from the electroweak scale. We refer to these LFV operators as “top operators”. 
Here, as an example, consider the exchange of a heavy SU(2) singlet leptoquark Sq with couplings XetSoeRt‘^ and 
X^cSofJ-Rc'^-, which (after Fierz rearrangement) generates the dimension -6 contact interaction 

^^^^(M 7 “^'fle)(c 7 “PRt) =-e^“‘^^( 7 l 7 “PRe)(c 7 “PRt) =. ( 1 ) 

Alternatively, we could define the operator coefficient as — 1 /A^, in which case e ~ ml/A? because 2-\/2 G'f — m^^ 
(we take mt = 173.3 GeV). We will quote low energy bounds on such interactions as limits on the dimensionless es. 
In the case of our leptoquark example, m^q ^ 1 TeV to satisfy current bounds on second generation leptoquarks 
from the LHC [22,23], and thus, for A ^ 1, one expects ^ 0.02. 

Notice that t and t —^ are mediated by different operators. Most of the bounds we quote 

will apply to _|_ |ge/ict| 2 ^ .^g study the LHC sensitivity assuming equal rates for t qe ~and 

t —>■ qe^ . 


2.1 Decay of the top 

In the Standard Model, the top decays almost always to 6 IF+, at a tree-level rate given by: 


T{t bW) = 




1 - 




l-b2- 


m'i 


~ 1.3 GeV 


In the presence of the operator of Eq. (1), the three-body decay rate is: 


( 2 ) 


r(t —>• + c) = \e 


fleet 1 2 GpTn^ 

' 1927r3 


(3) 


so the branching ratio, allowing for all the operators listed in the Appendix A, and neglecting fermion masses other 
than the top (to remove interferences), is 
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(4) 


where q G {u,c}, and we approximated yt\Vtb\ — 1- This is small ( 45 ^ — 2.8 x 10 due to three-body phase 
space. For the leptoquark example discussed above, BR{t —ce+/r“) ~ 10“® for = 0.02. 

The phase space distribution of the and q depends on the Lorentz and spinor structure of the contact 

interaction, and could affect the efficiency of experimental searches for this decay. The squared matrix-elements for 
the individual contact interactions have the form |A4p oc x(l — x), where x = and is the invariant 

mass-squared of a pair of final state fermions a and b. Our study will not take this into account, since we found, in 
some explicit examples, only a small relative effect on the selection efficiency (of the order of 5%). 


2.2 Bounds from flavour physics and HERA 

Low energy constraints on four fermion operators involving two leptons and two quarks have been estimated and 
compiled for many operators taken one at a time [24-26], and carefully studied for selected flavour combinations 
(see e.g. [27], or global fits to bsiIJI operators [28]). However, even in the more recent compilations [25,26], bounds 
on LFV operators involving a single top are not quoted. In the Appendices B and C, we estimate bounds on such 
operators from their possible contributions, inside a loop, to rare y, B and K decays. In Appendix D, we estimate 
bounds from single top searches at HERA. Here, we summarise the resulting bounds, and list in Tables 1 and 2 the 
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best limits on the coefficients of the various operators. We will find that only the coefficients of some operators are 
stringently constrained, while others could mediate LFV top decays within the sensitivity of the LHC. 

The current upper limit BR{^ —>■ ey) < 5.7 x 10“^^ [2] severely restricts e ^ flavour change. For our top 
operators to contribute, the quark lines must be closed, which requires at least two loops and a CKM factor, see 
the second diagram of Figure 5. Nonetheless, in the case of scalar or tensor operators involving Ir this diagram 
can overcontribute to /r —>■ ey by several orders of magnitude, because the lepton chirality flip is provided by the 
operator (rather than m^), so the diagram is enhanced by a factor We make order-of-inagnitude estimates 

in Appendix B, and quote the resulting bounds in Table 2. 

Exchanging a W between the t and q quark legs of the top operator will generate an operator with down-type 
external quark legs, see the left diagram of Figure 1. The coefficient of this light quark operator will be suppressed 
by a loop, CKM factors, and various masses. Numerical values for these suppression factors are given in Table 9 of 
the Appendix; however, their approximate magnitude is simple to estimate. If the top is singlet (tn), then the loop 
is finite; in the case of V zt A interactions, this is because mass insertions are required on both internal quark lines 
to flip chirality. In the case of scalar operators, one internal quark mass for chirality flip is still required, then terms 
linear in the loop momentum vanish, so the diagram is also proportional to an external quark mass. For scalar 
operators involving t r (which require an mq insertion to connect the W to the q line), the best limit can arise from 
exchanging a W between the t and a charged lepton leg, which generates a charged-current operator as represented 
in the second diagram of Figure 1. In the Appendix are also given current bounds on the coefficients of the various 
light quark operators that the top operators can induce. Comparing these bounds to the induced coefficients, gives 
the limits of Tables 1 and 2 that are labelled with (e)s. 




Figure 1: Diagrams for generating an LFV operator with light external quark legs, by dressing a top operator with 
a W loop. To reduce index confusion, down-type quarks have greek indices. 

HERA collided protons with positrons (or electrons) at a centre of mass energy of 319 GeV, and searched for 
single tops in the final state. The HI collaboration had a few events with energetic isolated leptons and missing 
energy, consistent with e^p — te^ + X followed by leptonic decay of the top [29,30]. However, this signal was not 
confirmed by the ZEUS experiment [31], and neither collaboration had a signal consistent with hadronic top decays. 
Both collaborations set bounds on a{e^p —>■ e^tV); we follow HI, since they had some events and a weaker bound: 

2 3 X10~® 

a^e'^p ^ + X) < 0.30 pb = -5- at 95%C'L. (5) 

mf 

Contact interactions of the form (eF/i)(uFt) and {'pTe){urt) could respectively mediate e~p —>■ p~tX and e~^p —>• 
fX^tX. As discussed in the Appendix D, the limit of Eq. (5) translates to e ^ 0.3 —>■ 1 for the various operators, as 
given in Tables 1 and 2. 

It can be seen from Tables 1 and 2, that rare decays give very weak bounds on some contact interactions of the 
form (eFT)(urt) and (rre)(?lFt). Such interactions might have contributed a signal at HERA via leptonic r decays, 
so we make some approximate estimates in Appendix D.4, and include the bounds in the tables with a 

2.3 Implications 

The current bounds on LFV branching ratios of the top can be obtained from Tables 1 and 2, and Eq. (4). In 
these tables, the bound on e^xYZ ^-ppears on the first row of each box, where the flavour indices ijqt are given in 
the left column, and the operator label XYZ is given in the first row (see Appendix A for operator definitions). In 
parentheses below the bound, is a clue to where it comes from: HERA means the single top searches at HERA that 
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Table 1: Constraints on the dimensionless coefficient of the four-fermion operator = 

2y/2GF {ei'yf^PxBj) {Uq-f^Pyt) ( see the operator list in Appendix A), where XY are given in the top row, and 
generation indices ijqt are given in the first column. The bounds are on the first line of each box, their origin 
on the second. They can arise from the HERA single-top search, or from the loop contribution to the operator 
involving lighter fermions whose coefficient is given in parentheses below the bound (see Appendix C.l for current 
experimental bounds on the lighter-fermion operator coefficients). The ~ bounds are discussed in Appendix D.4. 
We expect that e < 1, so in boldface are the bounds that impose an upper limit smaller than 1. 
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Table 2: Constraints on the dimensionless coefficient eb?*, of the scalar and tensor four-fermion interactions. See 
Appendix A for operator definitions corresponding to the subcript of e. The generation indices ijqt are given in the 
first column. Beneath each bound is given its origin in parentheses; £^s±p x limits of Table 7, and are 

from Table 8. See the caption of Table 1 for additional details. 


are discussed in Appendix D, and the (e)s mean the bound comes from dressing the top operator with a W loop. For 
instance, the bound 0.0037 in the second row and third column of Table 1, arises by exchanging a W between the 
t and u legs of the top operator, which gives the operator {ejfPLfJ-){djpPLs). This would mediate the unobserved 
decay 77° —>■ e/r, so its coefficient is bounded above, as indicated in Table 6. (The bounds on lighter-quark operators 
relevant to constraining top operators are given in Tables 6-8, and Table 9 gives the loop suppression factors with 
which the top operators generate the lighter-quark operators). Translated back to the top operator, the upper limit 
on BR{K° eqi) gives the quoted limit on the top operator coefficient. 

In this paper, we are interested in top decays to the bounds on which are given in the second and third 
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rows of the tables. For many^ of the operators involving the doublet component of the top {tL\ recall that the last 
index in the operator label is the top chirality), the rare decay bounds are restrictive, implying that these operators 
could only induce BR[t —> < 10“®. Scalar and tensor operators involving e, qi and t/j would overcontribute 

to ^ > ey. However, there remain operators which are weakly or not constrained, allowing a branching ratio ^ 10“^. 
It is therefore interesting to explore the sensitivity of the LHC to t —>• decays. 

Finally, it is interesting to consider how large the e coefficient of the top operators can be. Some of the upper 
bounds quoted in Tables 1 and 2 are ^ 1, and should not be interpreted as relevant constraints Indeed, the 
width of the top is given by DO [32] as 2.0±0.5 GeV (the theoretical decay rate to hW is 1.3 GeV), which constrains 
e^XYZ < 10-20. Furthermore, phenomenological prejudice and the leptoquark example of Eq. (1), suggest that e < 1, 
because the three-body decay should be mediated by sufficiently heavy (m > rrit) particles, with sufficiently small 
couplings to have not yet been detected. We therefore quote in bold face the “relevant” bounds that impose e < 1. 

3 t — e^jjL^q at the 8 TeV LHC 

In this Section, we estimate the sensitivity of current LHC data to the LFV top decays t where q = u,c. 

We consider strong production of a tt pair, because this is the most abundant source of tops at the LHC, followed 
by the leptonic decay on one top, and the LFV decay of the other. This is illustrated in Figure 2, and gives a final 
state containing three isolated muons or electrons which has small Standard Model backgrounds. 

3.1 Simulation setup 

This study is performed for proton-proton collisions at the LHC, with center of mass energy of 8 TeV and an 
integrated luminosity of 20 fb“^, corresponding to the LHC Runl. The details about the signal and background 
generation are given in Section 3.2. The detector simulation is carried out by Delphes [35] using a CMS setup 
parametrization. 

Delphes uses particle-flow-like reconstruction. The relative isolation of leptons is calculated from the total pr of 
the particles inside a cone of AR around the lepton direction (Ai? = 0.3 for electrons and 0.4 for muons), divided by 
the pt of the lepton. Jets are clustered using the fastjet package [36] with the Anti-kt [37] algorithm with distance 
parameter i? = 0.5. The b-tagging performances are tuned on the typical efficiency and fake rate obtained in CMS. 

For this study, no additional interactions in the same or neighbouring bunch crossing (pileup) are simulated. 

3.2 Signal and SM backgrounds generation 

The signal is generated with PYTHIA 8.205 [34] using tune 4C. Top quarks are pair produced, then one top is 
forced to decay to charm, , and e^, with equal probability between e~ and p~ e"*". The decay products are 
distributed according to the available phase space. 100k events have been generated both for LFV top and anti-top 
decays. 

The backgrounds for this search, listed in Table 3, are processes that can give rise to three isolated leptons and 
at least 2 jets in the final state. Most of them are related to the production of real isolated leptons, e.g. from a 
top pair or vector bosons in the final state. In the table, are also shown the details about the number of generated 
events and production cross section for 8 TeV and 13 TeV proton-proton collisions. The number of generated events 
refers to the generation at 8 TeV. The ti cross section is calculated with the Top-|--I-2.0 program to next-to-next-to- 
leading order in perturbative QCD, including soft-gluon resummation to next-to-next-to-leading-log order (see [38] 
and references therein), and assuming a top-quark mass of m* = 173.3 GeV. When an explicit calculation was not 
available, the cross sections have been calculated with the MCFM package [39], version 7.0. The kinematic cuts 
used for the calculation are also shown in the table. 

The leading order (LO) matrix element generator, MADGRAPH 5 [42], with CTEQ6 parton distribution func¬ 
tions, is used to generate top pair production, and associated production of a top pair and a vector boson (ttVF, 
ttZ). MADGRAPH, interfaced with tauola for r decays, is used to generate vector-vector production (VFIF, WZ 
and ZZ) and the contribution of weak processes giving rise to final states with one top quark, one h quark and a 
Z boson (decaying to leptons). For the vector-vector production, we only consider final states with at least two 

^The exception is the S — P operator involving whose loop suppression factor would involve two light quark masses. 

^They are given so that in the future, if the experimental bounds improve, the limits can be obtained by simply rescaling the number 
in the tables. For instance, if the upper bound on B t£X decays were to improve by two orders of magnitude, the limit on some es 
would be divided by 10, and become marginally relevant. 

®The final states where the W decays to t~Ut are not directly targeted by this search. The fact that such processes, followed by 
leptonic r decays, can pass our selection is taken into account in the signal efficiency, as explained in the following. 
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Figure 2: Feynman diagram for the considered signal where i 


e or /r (the conjugate diagram is also considered). 


Process 

Events 

a [pb] (8 TeV) 

cr[pb] (13 TeV) 

Source 

tt ( 212 i/ 2 b) 

12M 

26.19 

86.26 

Top-h-f2.0 (NNLO) [38] 

WW-fjets (212^, mu > lOGeVlc^) 

IM 

5.84 

11.57 

MGFM [40] 

ZZ-l-jets (212q, mu > lOGeV/c^) 

2.5M 

2.71 

5.35 

MGFM [40] 

ZZ-fjets (212 j^, mu > lOGeV/c^) 

IM 

0.774 

1.53 

MGFM [40] 

ZZ-fjets (41, mu > lOGeV/c^) 

2.5M 

0.390 

0.738 

MGFM [40] 

WZ-fjets (212q, mu > lOGeV/c^) 

IM 

2.37 

4.60 

MGFM [40] 

WZ-fjets {3lv, mu > lOGeV/c^) 

IM 

1.15 

2.23 

MGFM [40] 

ttW -fjets 

IM 

0.212 

0.612 

MGFM [40] 

ttZ-l-jets 

IM 

0.192 

0.798 

aMC@NLO [41] 

tbZ-fjets 

2M 

0.014 

0.047 

aMC@NLO [41] 


Table 3: Number of events generated, and cross sections at NLO (except for tt+jets), for each background category. 
Here I = e, fi, r. 


real charged leptons. This means that for the WW system, the considered final states are 2 charged leptons and 
2 neutrinos; for WZ, they are 3 charged leptons and one neutrino or 2 charged leptons and 2 quarks; and for ZZ^ 
they are 2 charged leptons and 2 neutrinos, 2 charged leptons and 2 quarks, or 4 charged leptons. In all cases, 
MADGRAPH accounts for the presence of up to two additional jets at matrix-element level, and the hadronisation 
is carried out by PYTHIA 8.205. The details about the SM backgrounds simulation and cross sections are shown 
in Table 3. 

3.3 Event selection 

The signal for this search is tt production, followed by the lepton flavor violating decay of one top (which will be 
denoted as LFV top in the following), and the leptonic decay of the other (standard top, in the following). The 
number of expected signal events is given by: 

Nsig = C ■ esiG ■ 2 • att ■ BRiW -)> liy) ■ BR{t -)> , (6) 

where I € {e,/r, r}, £ the integrated luminosity, esiG the selection efficiency on the signal, and att = 246.7 pb 
(see [38] and references therein). 

The considered signature is 3 isolated leptons (with one pair of opposite sign and opposite flavor from the LFV 
decay), 2 jets (one of which is a 6 -jet), and missing transverse energy. For the event selection, we consider only 
muons of pr > 20 GeV and jyyj < 2.4, electrons of pr > 20 GeV and jyyj < 2.5, and jets of pr > 30 GeV and 
jr^l < 2.4. These criteria are comparable to those used in real analyses by the GMS or ATLAS collaborations. A 
muon is considered isolated if its relative isolation value is less that 0 . 12 , and an electron is considered isolated if 
its relative isolation value is less that 0.1. We select events containing: 

1. exactly 3 isolated charged leptons (electrons or muons), two of which must be of opposite sign and opposite 
flavor. 
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Process 

no selection 

step 1 

step 2 

step 3 

step 4 

step 5 

step 6 

Signal 

202.57 

32.98 ± 0.17 

22.66 ± 0.14 

9.12 ± 0.09 

8.20 ± 0.09 

7.50 ± 0.09 

3.75 ± 0.06 

tt 

542806 

14.78 ± 0.81 

10.51 ± 0.69 

4.36 ± 0.44 

4.36 ± 0.44 

3.55 ± 0.40 

0.63 ± 0.17 

WW±jets 

116760 

0.93 ± 0.33 

0.35 ± 0.20 

< 0.35 

< 0.35 

< 0.35 

< 0.35 

ZZ-|-jGts 

72900 

353.74 ± 0.95 

82.50 ± 0.47 

3.74 ± 0.10 

1.60 ± 0.07 

0.25 ± 0.03 

0.03 ± 0.01 

WZ-l-jets 

63360 

852.21 ± 4.04 

182.96 ± 1.90 

8.70 ± 0.42 

7.62 ± 0.39 

0.74 ± 0.12 

0.04 ± 0.03 

tiw 

4240 

9.36 ± 0.24 

7.67 ± 0.22 

3.59 ± 0.15 

3.45 ± 0.14 

3.10 ± 0.14 

0.27 ± 0.04 

tiZ 

3840 

17.25 ± 0.33 

16.44 ± 0.32 

7.72 ± 0.22 

7.16 ± 0.21 

1.85 ± 0.11 

0.22 ± 0.04 

tbZ 

282 

5.75 ± 0.03 

3.59 ± 0.02 

1.51 ± 0.01 

1.37 ± 0.01 

0.13 ± 0.01 

0.01 ± 0.01 


Table 4: Number of expected events for a luminosity of 20 fb~^, at various steps of the selection, for the signal 
process normalized to a branching ratio BR(t —?► qe/r) = 6.3 • 10“®, and the various backgrounds considered in this 
study normalized to their NLO cross sections. All uncertainties are statistical only. The considered backgrounds 
are the same as in Table 3, grouped in wider categories. In particular, the numbers in the ”no selection” column 
are relative only to the final states detailed in Table 3. The steps in the selection are, step 1: 3 leptons with two 
opposite sign opposite flavour, step 2: at least 2 jets, step 3: exactly one b-tag, step 4: missing Et greater than 20 
GeV, step 5: Z boson veto, step 6 : invariant masses cuts (see Section 3.3 and Figure 3). Uncertainties are statistical 
only. 

2. Events are requested to contain at least 2 jets, and 

3. exactly one b-tagged jet, and 

4. the missing transverse energy has to be higher than 20 GeV. 

5. In order to exclude events where two of the isolated leptons come from a real Z boson, we reject events 
containing any pair of opposite sign isolated muons or electrons with invariant mass between 78 and 102 
GeV/c^. This cut is particularly helpful in rejecting background arising from tt associated production with a 

Z. 

6 . The charged lepton that does not belong to the pair of opposite sign and opposite flavor leptons is assigned 
to the standard top in the event, and assumed to come from the W decay (bachelor lepton). Following 
a common procedure in reconstruction of tt semi-leptonic events, the x and y components of the missing 
transverse energy are taken as a measurement of the neutrino Px and py, and the longitudinal component 
of the neutrino momentum is calculated imposing that the invariant mass of the system composed of the 
bachelor lepton and the neutrino must equal the mass of the W boson. The bachelor lepton and the neutrino 
4-momenta are then combined with that of the b-tagged jet, to build a candidate standard top. When more 
choices of the bachelor lepton are possible (there can be up to 2 possible pairs of opposite sign opposite flavor 
charged leptons in one event), all are considered and the one giving the best standard top mass is chosen. We 
reject events in which the invariant mass of the standard top candidate is more than 45 GeV away from the 
nominal top mass. After the choice of the bachelor lepton, there is only one possible pair of opposite sign and 
opposite flavor leptons in each event. This is combined with all good (non b-tagged) jets present in the event 
to build a list of candidates for the LEV top. 

Events are requested to have at least one combination giving a LEV top mass within 25 GeV of the nominal 
value. 

The efficiency of the final selection on signal events is 

es/G = (1.85 ±0.03)% , (7) 

where the uncertainty is statistical only. The signal efficiency is calculated on tt events where one top decays through 
t —>• and the other one decays to a 6 quark and a W, which subsequently decays to a charged lepton (e, p 

or r) and a neutrino, and is defined as the fraction of such events passing the selection criteria. 

The number of expected events, for the signal and for each background category, on 20 fb~^ of proton-proton 
data at 8 TeV is shown in Table 4, for different subsequent selection requirements. 
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3.4 Results and expected limits on the branching ratio 

The selection and its efficiency, on signal and background, are discussed in Section 3.3, and summarised in Table 4. 
Assuming a branching ratio of BR{t qn^e^) = 6.3-10“^ for the signal, an uncertainty of 2.5% on the luminosity, 
and 20 fb“^ of data, we would expect Nsig = 3.75 ± 0.06 signal events, to compare to the Nbkg = 1-20 ± 0.18 
expected events from known backgrounds. 

In Figure 3, we show the invariant mass of the LFV top candidate (left) and the standard top candidate (right), 
in events passing all the cuts except those on the masses themselves. 
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Figure 3: Invariant mass of the LFV top candidates (left) and standard top candidate (right) in events passing 
all the selection, apart from the cut on the masses themselves. The different background contributions are shown 
in filled histograms and stacked, normalised to the number of events expected in 20 fb“^. For comparison, the 
distribution for signal events is also shown as a dashed line, normalised to the number of events expected in 20 
fb“^, for an example signal branching ratio BR(t —>■ = 6.3 • 10“®, equal to the expected limit extracted in 

Section 3.4. 

In order to evaluate the sensitivity of this search for t —?► we calculate the expected upper limit that 

could be set, in the case of absence of the signal. The calculation is based on the number of expected background 
events surviving the final selection, in 20 fb“^ of 8 TeV LHC data, so the result can be interpreted as the possible 
upper limit the CMS or ATLAS collaborations (alone) could expect to set with RunI data, if the signal were not 
there. 

A 95% confidence level (CL) upper limit on the branching fraction of t —>■ qfi^e^ is calculated using the modified 
frequentist approach (CL^ method [43,44]), as it is implemented in the RooStats framework [45]. Based on the 
number of expected background events, summarised in Table 4, and on Equation (6), the obtained limit is: 

BR{t qfi^e^) < 6.3 • 10"^ @95%CL. 

Alternative techniques for limit calculations, as implemented in RooStats, have been tried, leading to compatible 
results. An eventual variation of 100% in the number of expected background events would lead, in the worst case, 
to an expected limit of BR{t q^^e^) < 7.4 • 10“® @ 95 % CL. 

As explained in appendix E, we emulate in our framework the published CMS search for BR{t Zq) [11]. 
The reason for this exercise is twofold: on one hand it allows to validate our procedure on simulated samples, by 
comparing with the CMS background expectations. On the other hand, it provides an estimate of the constraint set 
on BR{t —>■ qfi^e^) from this previous analysis, which is found to be BR(t —>■ qfi^e^) < 3.7 • 10“^, on the verge of 
probing LFV top decays mediated by a four-fermion operator. As proven in the present study, a dedicated analysis 
would set a limit 50 times stronger (of the order of BRit —^ qfi^e^) < 6.3 • 10“^), showing that the existing LHC 
data from Runl can still be used to obtain interesting constraints on lepton flavor violation. 
























Figure 4: Single top production via the LFV contact interaction, which produces a final state similar to the 
considered signal. 


4 Discussion 

4.1 Perspectives at 13 TeV and 14 TeV 

To estimate the reach of the described search at a center-of-mass energy of 13 TeV, we extrapolate the 8 TeV 
results, rather than performing a full simulation of signal and background processes at 13 TeV. The increase of 
the production cross sections for SM processes, from 8 TeV to 13 TeV (see Table 3), is taken into account. The 
selection requirements and efficiencies are kept the same as for the 8 TeV analysis. For the signal, we have checked 
on simulated events that the efficiencies at 8 TeV and 13 TeV are consistent within 5%. 

The sensitivity is estimated by calculating the expected upper limits on BR{t — qn^e^), in the absence of 
signal, for two scenarios: the case of 20 fb“^ and 100 fb“^ proton-proton data collected by the LHC at 13 TeV 
center-of-mass energy. We also extrapolate the sensitivity to the case of 3000 fb^ of integrated luminosity at 14 
TeV. In this last case, we simply rescale signal and background rates from 13 TeV to 14 TeV, and use the square 
root of the number of expected background events as an estimate of their uncertainty. The obtained values are 
summarized in Table 5. The upper limits presented here are derived using statistical uncertainties only, so don’t 
take into account the possibility for such analyses to become systematically dominated in the future. In order 
to have an accurate evaluation of the systematics evolution, a deeper study from the LHC experiments would be 
needed. On the other hand, for an analysis on 13 TeV or 14 TeV data, the selection would have to be re-optimised, 
possibly leading to an increase in sensitivity. 



8 TeV (20 fb-") 

13 TeV (20 fb-") 

13 TeV (100 fb-") 

14 TeV (3000 fb"") 

BR{t -A qpL^e^) 

< 6.3-10-'^ 

< 2.9-10-^ 

< 1 . 2 - 10 -'^ 

<2 - 10 -“ 


Table 5: Expected upper limits on BR{t —>■ in the hypothesis of the absence of signal, for 8 TeV, 13 TeV (in 

two scenarios: the case of 20 fb~^ and 100 fb“^ collected luminosity) and 14 TeV for 3000 fb“^ collected luminosity. 


4.2 Single Top 

In addition to mediating LFV top decays, the top operators listed in Appendix A could lead to single top production 
with an pair, as illustrated in Figure 4. The objects in the final state would be the same as for the tt process 

we studied : three leptons, missing energy, and > 2 jets, of which one is a b. We estimate^ that at 8 TeV, the cross 
section for pp —>• —>■ e^pL^,Jvb is of similar order to the cross section for pp ^ it ^ e^p^q,Jvb, for operators 

involving a u quark and slightly less for a c quark. 

^The partonic diagrams for u + g ^ e~fB't are those oi b g ^ W + t, with W replaced by par e~/r+, and 6 by « or c. So we 
estimate (7{pp e~p^t) by comparing to cr{pp —>■ Wt) ~ 24 pb ~ 0.1 X cr(pp —>■ tt). If we neglect spin correlations between the e~ 
and ut, then for V ± A top operators, \M(ug —>■ e“/r+t)p can be factorised as oc ^\A4(ug —>■ Z't)\'^ X \M(Z' —>■ e~p'^)\'^ where the Z' 
has coupling 1/A to fermions, and the 1/3 is the spin average. Then one has \M{Z' e~p'^)\'^ = 2Q^j where = (pe -t-p^i)'^, 

and \M{ug —>■ Z't)\'^ = 2\M{bg —>• Wt)\^/(h?g'^). The three body phase space integral can be written as J dHyyt X dlle^ X dQ'^ 1(211), 
where dllyvt integrates the phase-space of W and t in the centre-of-mass frame, J dUe^ = l/(87r), and we take the upper limit of the 
integral as assuming thhis is a reasonable value for leptons with pT > 20 GeV. Finally, at mt/(8 TeV) ^ x ^ 2mt/8 TeV, the 
ratio of tt to 6 pdfs is of order 6-10, so we estimate that cF(pp —^ ep'^t) ~ cr(pp —^ ti) X BR(t —> qp^e^). 
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We neglect this process for two reasons. First, the contact interaction approximation (that the |four-momentump 
in the process <C A^) is more difficult to justify than in top decay, because the energy scale in the process can be 
^ rrit. Secondly, finding such events in the backgrounds could be more challenging because the cannot be 

required to participate in the reconstruction of a top (step 6 of Section 3.3). 

We envisage that it makes sense to neglect the LFV single top process in a first search for LFV top decays. This 
is conservative, because LFV single top production could contribute events that pass our selection.In the absence 
of a signal, such a search could sufficiently constrain the contact interaction scale A, to justify including the single 
top process in subsequent analyses. 


5 Summary 


The aim of this paper was to explore the LHC sensitivity to the decay t —>■ which is lepton and quark 

flavour-changing, but baryon and lepton number conserving. We parametrise this decay as occuring via a contact 
interaction, and list a complete set of SUiZ) x 17(1) invariant dimension-6 operators in Appendix A. We parametrise 
the coefficient of these interactions, which we refer to as “top operators”, as e2-\/2G_F or equivalently 1/A^, with 



Model-building prejudice (see Section 2.3), suggests that e ^ 1. The top branching ratio is then 


BR{t^qi+i-) ~ 1.3 


487r2 


1 V±A 

i S±P 


These contact interactions are currently constrained from their contribution in loops to rare decays, and from 
single top searches at HERA. These bounds are discussed in the Appendices, and summarised in section 2. For 
interactions involving e and /r, rare decay bounds impose e ^ 0.01 for some operators, but others can have e ^ 1. 

In Section 3, we evaluate the sensitivity reach of a dedicated search for lepton flavor violation in top decays, 
at the 8 TeV LHC. The search targets tt events, where one top decays to an up-type quark (u or c) and a pair of 
leptons of opposite sign and opposite flavor, and the other one decays to a 6 quark and a VF, which subsequently 
decays to a charged lepton and a neutrino. This is illustrated in Figure 2. 

The relevant signal and SM background processes are simulated for LHC Runl-like conditions: proton-proton 
collisions at 8 TeV center of mass energy, for an integrated luminosity of about 20 fb“^. The detector simulation 
is based on Delphes, with parameters tuned on the CMS detector reconstruction and performances, but does not 
include pileup. The analysis setup is validated by emulating an existing CMS search for rare top decays to Zq in 
tt events, showing reasonable results. 

We find that a dedicated search by a single experiment using 20 fb“^ of 8 TeV data could be sensitive to 


BR{t eqi+ jet) ^ 6.3 x 10 ® , 


and extrapolate that a sensitivity of ~ 1.2 x 10“® (^ 2 x 10“®) could be reached with 100 fb“^ at 13 TeV (3000 
fb“^ at 14 TeV). From Equation (4), we see that the 100 fb“^ data could impose |e| < 0.06 for V ± A operators, 
and |e| < 0.1 for the S' ± P and LQ operators. This analysis shows that the existing LHC data from Runl can 
still be used to obtain interesting constraints on lepton flavor violation. Although this is understandably not the 
priority focus in the times of the Run2 startup, let’s not to leave unchecked this possible path to New Physics. 


Acknowledgements 

SD thanks P Gambino and Professor Y. Kuno for useful conversations, and the labex OCEVU and the theoretical 
physics group at Montpellier University-H for hospitality during part of this work. The work of M.L.M. is performed 
in the framework of the ERC grant 291377 “LHCtheory: Theoretical predictions and analyses of LHC physics: 
advancing the precision frontier”. 


A Appendix: Operators 

Consider the S'-matrix element mediating t —>■ q£^ij, where £k € {e,/r, r}, and suppose it is induced by local 
operators with momentum-independent coefficients, such that they can be added to the Standard Model Lagrangian. 
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These operators should respect the SU{3) x SU{2) x 17(1) gauge symmetry of the Standard Model. However, since 
the New Physics scale that we explore is not much larger than the electroweak scale mt, we should include dimension- 
8 operators constructed from two Higgs helds and four fermions, or two gradients and four fermions. Instead, we 
choose to work with SU{3) x {7(l)-invariant, but not S'?7(2)-invariant, operators of dimension 6. This is because a 
dimension-8 S'{7(2)-invariant operator, can be a dimension-6 SU{3) x 17(l)-invariant operator: 

where Qi is a quark doublet of *th generation, Ej is a lepton singlet, = (i7+,i7o), = iT2H* and SU(2) 

contraction is in square brackets. 

This choice of S't/(2)-non-invariant operators means that the coefficients of (ejaPxPLb) and {e'yaPxfJ-)x 
(c 7 “Pit) are taken as independent, and in particular, bounds on the first do not apply directly to the second. 
However, dressing the top operator with a W loop will generate the light quark vertex, which gives unavoidable 
bounds that are estimated in Appendix C. (Treating the W as dynamical while the New Physics is a contact 
interaction is unlikely to be a good approximation, but is the only way we can estimate whether the top operators 
are in tension with other observables.) In the case of /x —>■ ey, two or three loops are required to transform the top 
operators into the dipole operator, so we need an SU(2)-invariant formulation of the top operators. These dimension 
six and eight operators are given in Appendix A.2. 

A.l SU{3) X t/(l) invariant “top operators” 

We are interested in contact interactions involving two charged leptons of different flavour, a top and a c or u quark. 
Such colour-singlet, electric charge-conserving, dimension-6 operators will be referred to as “top operators”. The 
V ± A operators are: 


= {ei'y°‘PLej){uq-faPLt) (9) 

= {ei'y°‘PLej){uqjaPRt) ( 10 ) 

= {ea^PRe,){uq-fM) (II) 

= (e,7“PKe,)(u,7aPxxt) (12) 


and a redundant list of scalar/tensor operators is: 


qS+P,R _ 

{eiPRej){uqPRt) 

(13) 

qS+P,L ^ 

{eiPLej){uqPLt) 

(14) 

qS-P.R _ 

{eiPLej){uqPRt) 

(15) 

qS-P.L ^ 

{e^PRej){uqPLt) 

(16) 

qT.R ^ 

{uqa^’' PRej){eia^„PRt) 

(17) 

0 Ti = 

{uqa^"’' PLej){eiafj,^PLt) 

(18) 

0LQ,R _ 

{uqPRe,){e.PRt) = -^O^+PLi + l0T.fl 

(19) 

^LQ^L _ 

{uqPLeq){e,PLt) = 

(20) 


where e, u are Dirac spinors, i,j are unequal lepton flavour indices, q £ {m, c}, and the last chiral superscript of the 
operators gives the chirality of the top. 

The scalar operators 0^*7A exchanged for the tensor operators , as shown in Equations (19) and 

(20). We will use the LQ and T operators interchangeably, becuase the 0^*3A g^pg jnore convenient in top decay, 
and the tensors at low energy. The tensors are used in the basis of [46]. Notice that S + P operators are defined to 
have the same chiral projector twice, whereas for S — P they are different. Since the last quark flavour label is fixed 
to S + P, S — P and LQ operators appear twice, for both chiralities and tn (were we using arbitrary quark 
flavour indices q,v, then Finally, the scalar operator contracting quark to lepton spinors is 

given the name LQ] there are no such operators of the form {uqP]iej)(ejPLUt), because this is Fierz-equivalent to 
qRL 

A similar list can be constructed for down-type quarks. These will be relevant, because the top-operators can 
generate the down-operators at one loop. To reduce index confusion, down-type quarks will have greek flavour 
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( 21 ) 


indices a,/3 S {d, s, 6}. There are also charge-current (CC) scalar operators 

= (l^PLVj){UqPLdp) , oCC.LQ ^ {uqPLUjmPLdp) 

which will be relevant in setting bounds. 

These operators are included in the Lagrangian with a dimensionful coefficient 

-^ = - 2 \/ 2 GFe . ( 22 ) 

The es have subscripts that identify the operator, and lepton-quark flavour superscripts: 2\/2GFilL would be the 
coefficient of PLej){uq"faPLt), and ^Iq x would be the coefficient of (uqPxej){eiPxt) (so in particular, the 
index order is always liqq, even for O^^). 

A.2 SU{3) X SU{2) X (7(1) operators 

Our V ± A operators can be identified as components of SU(2) invariant dimension 6 operators: 


qLL ^ (e7“PF/r)(lIg7„Pit) C (7^ = (Le7“iF)(Q,7cQt) (23) 

= [(e7“PL/i) -f (T7e7“PLt'^)][(Mg7ar“PLt) -k (dqjaT°-PLb)] 

qLR ^ {erPLt^){uqlc.PRt)cOLu = (Lel^L^){Uq^^Ut) (24) 

= ie-t‘^PRfi)(uqjc.PLt)cOEQ = (Eej‘^E^){QgjM (25) 

qRR = (ej‘^PR^i){uqj^PRt) = OEU = (Eel‘^E^)(UqjM) (26) 


where L, Q are doublets and E, U are singlets t is a third generation quark index, q a first or second generation 
quark index, and we explicitly choose lepton indices i = e and j = /r. In the SU(2) invariant list, there is an 
additional operator at dimension six, which generates Oll in combination with charged current interactions: 

= (lG“T“L,)(Q, 7 „T“Qi) 

= 2{Vi^°^PLej)(dq-faPLt) + 2{eiJ°‘PLRj){Uq^aPLb) 

+ [{Vi'J°‘PLRj) - {e^'y°'PLej)][{uqjaPLt) - (dqjaPLb)] ■ (27) 

Alternatively, at dimension eight, SU(2) invariant operators can be constructed to contain only an up-type quark 
current, by contracting quark and Higgs doublets. An example is given in Equation ( 8 ). 

Then there are two S + P operators: 

qS+p,r _ {ePR^){uqPRt) C OL 1 E 2 Q 2 U 3 
qS+PA ^ ^ 

which are components of the same SU(2) invariant operator: 

Olequ = {Li Ej)eAB{Qq Ut) (28) 

= -{piej)(dqUt) -I- {e%PRej){uqPRUt) 

and there are the LQ or tensor operators, which are components of the dimension six, SU(2)-invariant operator 

Oluqe = {Li Ut)eAB{Qq Ej) = --Olequ + -^{Li a'^‘'Ej)eAB{QqCrfj., 2 Ut) , (29) 

and finally there are two S — P operators, which can be identified with an SU(2)-invariant operator of dimension 
eight: 

V^O^-P'R = V^{ePLfJ.){UqPRt) (30) 

= {E4H^L2]){[QqH^]U,) 

V^OS-P’L = {[LeH]Ei,){Uq[HPtT 2 Qt]) (31) 

where, for instance, [iJ^L 2 ] = vfiL- Notice that our S + P and tensor operators can also be obtained at dimension 
8 without associated charged current interactions: 

{L,HE,){QgH^Ut) = v^e,PRej){uqPRUt) (32) 
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B Appendix: top operators in /i ^ 57 

In this Appendix, we wish to close the quark lines of a top operator, and attach a photon (and a Higgs vev) to the 
resulting diagram, such that it contributes to eq. These two and three loop estimates rely on using SU(2)- 
invariant operators (of dimension six or eight), because using the di, and components of doublets will allow to 
avoid GIM suppression. Then we estimate the contribution to /r —?► ey by power counting selected diagrams. 


B.l Parametrising n ^ e'y 

The decay /x —>■ ey is mediated by the dipole operator, which can be added to the Standard Model Lagrangian as [1] 


£ 


meg 


AGf 

■ vf 

AGf 

'TT 


rri/j, {A[i^ji(^°‘^eLFap + CFiFap) 

Vfi {Aii]I^a°‘^[F{^Le]Fai3 + AL[L^F[]a°'^enFap) 


(33) 


where, in the second formulation, is the muon Yukawa coupling such that = m^, and the fermion part has 
been written in SU(2)-invariant form, to emphasize that the dipole operator is of dimension six, and has a Higgs 
leg. This operator gives a branching ratio: 


BR{y ey) = SSAtt^HArI'^ + \Al\'^) < 5.7 x 


(34) 


where we quote the upper bound from [2]. If \Ar\ = \Af\, this implies \Ax\ < 8.6 x 10“®. However, in the case of 
our diagrams, there will be a sum over three colours at the amplitude-squared level, so its convenient to divide by 
■\/3, and impose 

< 5.0 X 10“®. (35) 

This is a small number. We can estimate the contribution of top operators to be suppressed by various factors: 


1. two loops, and in some cases three. Two loops are neccessary to close the quark legs, because they are of 
different flavour, so must interact with a W. 

2. CKM, because the quark legs are of different generation. This is ~ .04 for cot flavour change, ~ .008 for 
u O t flavour change. 

3. GIM: if the W is exchanged between two up-type quark lines (see Figure 5), then the diagram is further 
suppressed by a GIM factor ^ rnl/m^ ~ 1/400 (for mf, = 4 GeV). This can be avoided by exchanging the 
W between the quark and lepton lines. 


4. mass insertions: may be required on the quark line to transform singlet into doublet quarks. However, the 
dipole operator of Equation (33) is defined with a Higgs leg attached to the muon line; so this can be an 
amplification factor if the Higgs can be attached to a heavier fermion. 

So for our SU{3) x ?7(1) invariant top operators with coefficient one could anticipate 


Ax 


egV2 

( 1671^)2 


e X GKM X GIM x m—insertations ~ 


10 ^ ex GIM X m — insertations 
2 X 10“® e X GIM x m — insertations 



(36) 


where, in some diagrams, the GIM factor is absent, and/or the mass insertion factor can be ^ 1. Gomparing to 
Equation (35) suggests that /x —> ey can provide relevant constraints on the es, unless there is more suppression 
than CKM and two loops. 


B.2 Estimated contributions to yu —)■ 07 

It is convenient to pretend that the New Physics scale is ^ v, which allows to study whether (the SU(2) invariant 
generalisations of) the top operators mix to the dipole in electroweak running down to v. The aim is to add 
electroweak boson loops to the top operator such that the quark lines can be closed, and the dipole operator 
of Equations (33) is obtained. Closing the quark lines requires changing the quark flavours. The lepton spinor 
contraction of the dipole is tensorial, so for all non-tensor top operators, it must be changed to this form by 
involving a lepton line in a loop. 
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1. Consider first the SU(2) invariant versions of the ViLA operators of Equations (23-26). These operators have 
no {e'yu)(u'yd) components, so two W vertices are required on the quark loop (see the first diagram of figure 5). 
This gives 2-loop and CKM suppression. For singlet quark currents, also mass insertions are required on the 
t and q lines, and the diagram will be suppressed by a GIM factor ~ m^/rn^ (these mass factors correspond 
to the crosses on the quark loop of the first diagram of Figure 5). To modify the spinor contraction between 
the leptons, a 7, Z could be exchanged between the quark loop and an external lepton line, which gives an 
additional suppression ~ e^/(167r^). So for singlet quark currents (the operators of Equations (24) and (26)): 


Ax 


eV/2 


€ X Vgb 


ml mg 
mly mt 


~ 10"^^e X 


.04 


(37) 


and for doublet quark currents, where the GIM and mg suppression can be avoided by emitting d-type quarks 
from the operator: 


Ax 


e^gV 2 

(167r^)3 


e X Vtg 


X 


.04 


(38) 


These loops gives no bound on ell, ^lr, C-rl or tRn- 




Figure 5: Left: example three-loop diagram by which a the neutral current vector operator of Equation (24) could 
induce /i —>■ eq. “x” is a mass insertion, the grey blob is the operator, and the photon could be attached to any 
quark line. On the right is an example two-loop diagram by which a tensor/scalar operator could induce /i —?► eq 
without GIM suppression and with mt/m^ enhancement. 


2. However, if the Oll operator is generated as a component of the triplet SU(2)-invariant operator of Equation 
(27), then there is a charged current contact interaction, which allows W exchange between the lepton and 
quark lines and avoids any GIM suppression. The diagram is similar to Figure 5 on the right, but with 
and Cl instead of and cr, and a mass insertion on the external /r rather than the t. An estimate gives 


Tlx 


egV 2 

(167r2)2 


e X CKM ~ I0"^e x 


CKAI 

.04 


(39) 


which imposes crl ^ 0.03 for g = c, or crr ft, 0.5 for q = u. Note that the B decay constraints, which are 
estimated below more rigorously than these ones, are tighter. 


3. Gonsider now the S' -I- T* or T operators which arise at dimension six, which can induce the second diagram of 
Figure 5. The grey blob represents the charge-changing component of the operator, involving tR and either 
d-L or SL- The contribution to ^ eq is two-loop and GKM suppressed, but enhanced by mt/m^, because 
the chirality flip for the leptons is in the contact interaction, so the Higgs leg of the dipole operator can attach 
to the mass insertion required to flip the chirality of the top. This gives 


Ax 


egV2 J 

' Vts 

(16x2)2 6'' 1 

^ Vtd 

^ lO-®! 



mt 

mu 


2 X 10-4 

6 X I0-® 


(40) 


which implies es+p,R, f-T,R ^ 10 ^ 

In the case where the third generation quark is a doublet, then there would be a bp and a charm or up mass 
insertion, so that 


A 


X 


egV 2 

( 16 ^ 2)2 

which would suggest r, ^t.r ^ 10“^. 


e X Vgb^ 
m,. 


10-6 
6 X 10-4° 


( 41 ) 
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4. Finally, consider an 5” ± P or T operator of dimension eight, so with two Higgs legs. If the Higgs are both 
neutral vevs, then the contact interaction involves up-type quarks, and gives diagrams whose contribution 
to /i —>■ ey is suppressed by three loops, CKM and GIM. Instead, focus on diagrams where the Higgs legs 
are charged, so the four-fermion interaction is charge-changing for the leptons and quarks, and the Higgs 
loop has to close (two external Higgs legs can be attached to the top line, to keep the operator at dimension 
eight, which gives negligible suppression because y* ^ I). Due to the Higgs loop, the estimate for Ax will be 
suppressed by ^ 10“^ with respect to the previous discussion of 5 -|- P and T at dimension six. Since we do 
not know whether our “top operators” correspond to SU(2) invariant operators of dimension six or eight, we 
conservatively retain the dimension eight bounds, which imply that for S' ± P and T operators involving a 
singlet tfi: 


e 


ct 

TM ’ 


^S±P,R 


^ 10-3 


ut 


ut 

^S±P,R 


^ 10-2 . 


(42) 


Similar estimates can be made for the contributions of top operators to t —>■ £ 7 . The current experimental upper 
bounds [53] are r(r —ey) < 1.8 x 10-'^r(T —>■ epp) and r(T —^ /ry) < 2.5 x 10 ^r(T —^ Upp)-, which imply, with 
the assumptions leading to Equation (35), that Ax{t —?► £y) ^ 3 x 10“®. For dimension eight, S ± P or tensor 
operators, this would give ^100 and ^10. 


C Appendix: top operators in meson decays 

Dressing a top operator with a W loop, where the W attaches to the top leg, generates a contact interaction among 
light fermions, where the t is replaced by 6,s or d. So the first step in obtaining bounds on the top operators, is to 
collect up-to-date bounds on the relevant two-quark-two lepton operators invoving d-type quarks. These are listed 
in Tables 6, 7 and 8, updated from [25]. In Section C.2 , we estimate the loop factors, which, when multiplied by 
the top-operator coefficients, give the coefficients of light-fermion operators that are constrained by Tables 6, 7 and 
8 . Finally, in Section C.3, we combine the loop factors and Tables to obtain the bounds on top operator coefficients 
given in Tables 1 and 2, and briefly discuss possible cancellations among the operators. 

C.l Updated rare decay bounds 


ijaP 


Constraint on 

Observable 

Experimental value 

efjids 

(cll — £lr) , i^RR — £rl) 

3.0 X 10-^ 

BRiK^ e^^) 

< 4.7 X 10-4^ 

efidb 

{epL — slr) 

f-LL + ^LR 

(epp — epp) 
epp + eRp 

3.0 X 10-4 

1.3 X 10-4 

BR{B^ ey) 

PP(P+-)'7r+e/i) 

BR(B'>^'K-ev) 

< 2.8 X 10-y 

< 1.3 X 10-3 

efisb 

Ell + CLP 

epp + f-RL 

1.0 X 10-4 

BR{B+^K+eil) 

, 9.1x10“® 

BR(B+^D°eve) 

^ 2.2x10-2 

erds 

(cLL — SLr) 

(epp — epp) 

4.1 X 10-4 

BR(T^eK) 

BR(t^vK) 

^ 2.6x10“® 

^ 7.0x10“® 

erdb 

{f-LL — ^Lr) 

(epp — epp) 

2.1 X 10-3 

BR{B^ er) 

< 2.8 X 10-3 

ersb 

ePL + CLP 

epp -b epp 

2.0 X 10-3 

BR(B+^K+fe) 

BR{B+^D°tu) 

^ 3.0x10“® 

^ 7.7x10“® 

firds 

(epL — ^lr) 

(epp — crl) 

4.3 X 10-4 

BRir^nK) 

BR(t^uK) 

^ 2.3x10“® 

^ 7.0x10“® 

^rdb 

(epL — £lr) , (epp — £rl) 

1.6 X 10-3 

BR{B^ fir) 

< 2.2 X 10-4 

^Tsb 

f-LL + ^LR 

epp + eRp 

3.1 X 10-3 

BR{B+^K^ffi) 
BR(B + ^D'>tu) 

^ 4.8x10“® 

^ 7.7x10“® 


Table 6: Constraints on the dimensionless coefficient of four-fermion interactions 

2 'J2 Gf {&il^PL,R^]) {da"fiJ.PL., Rdp). Pl,r can be Pl or Pr. The generation indices ijaf3 are given in the 
first column, and the best bound in column three. It arises from the observable of column 4, and the experimental 
value we used is given in column 5. In column 2, is given the linear combination of es to which the bound applies. 
All bounds apply under permutation of the lepton and/or quark indices. 


The bounds from the decay rate of a pseudoscalar meson M, are obtained using the usual formula: 


F (M ^ PP) = 


kGi 


‘‘M 


P^ [(|es±p,fl - esiP.il^) {m\i “ '<^1 “ ”^|) + 4eLeRmim2] -I- 
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ijafd 

Constraint on 

Observable 

Experimental value 

peds 

9.0 X 10-9 

BR{Kl fie) 

< 4.7 X 10-42 

epbd 

8.1 X 10-9 

BR{B'> eg) 

< 9.2 X 10-3 

ephs 

2.9 X 10-3 

BR{B'^ eg) 

< 1.1 X 10-9 

erds 

2.9 X 10-^ 

BR{t K'^ge) 

< 2.6 X 10-3 

erdb 

5.6 X 10-4 

BR{B° er) 

< 2.8 X 10-3 

ersb 




pLTbd 

5.0 X 10-4 

BRiB'^ fir) 

< 2.2 X 10-3 

jirds 

2.6 X 10-4 

BRir K^ge) 

< 2.3 X 10-3 

piTsb 





Table 7: Constraints on es±p.x, for the ijaf3 index combination given in the first column, where the es are coefficients 
of the operators Os±p,x = {eiPyej) {daPxdp), in a basis containing also tensor operators. The second column is 
the constraint, which arises from the observable given in column 3. The experimental value used is the last column. 
These bounds are also valid under lepton and/or quark index permutation. 


e-xR - [{m\i -mj - + to|) + 4TO-m|] \ , 


(43) 


where scalar-vector interference was neglected. In current algebra, the quark currents are: 


Ap^ = \mriWM) = 


5 ^/nl- 5 iA/r\ Imppm PPM 

P = x{moi qn\M) = —-— 

2 2 


(44) 


where o, n are flavour indices, and k is the magnitude of the final state three-momenta in the centre-of-mass frame: 

1 


= 


Am ^ 


(jPm - {p^M - (p^i - 


(45) 


The decay —>■ K~^e~£'^ can be used to constrain with a suitable approximation for the hadronic matrix 
element. For this, we follow the Babar exclusive determination of V/h [47] from —>■ where the differential 

decay rate is parametrised with a a kinematic term depending on = (pb—PdY ^ form factor, that is function 
of the inner product w of i? and D 4-velocities: 


duj 


G%\Vcb\'^m% 
8 X 487r3 


{niB + TOn)^(l + 0 {mjjjrp%) 




where ... includes the terms. To extract a bound on e from > K'^e r’*', we suppose the form factors are 
similar, and impose 

r{B+^ K+e-T+) {tpb+ppkY 

r(i?+ ^ DH+u) ~ IKbP {rPB + rnoY ’ ^ ’ 

which is slightly weaker than the bound in [25], for the case £ = fi. 

For S ± P operators, we included r — Ks£ bounds, that were not in [25], but used in [48]. 


C.2 Estimating the loops 

Consider first the V ± A operators involving the singlet with coefficient 2y/2GFeYR^ g S {c,u}. Dressing 

with a W loop as in Figure 1 gives the operators = {£i'yaPY£j) {da"t°‘PLdp) with coefficients 


ijaP 

^YL 


g^pptPPqVtfiV*c. , ( 


‘2‘OLem 

3 


VtpV*^ 


PPq 

PPt 


Jjqt 

^YR 


(47) 


where (/?, a) G {(6, s), (6, d), (s, d)}, g G {c,u}. The log is ~ 1.54 ~ tt/2. The numerical factor arising from the 
loop depends on the light external quark indices /3, a, as well as the quark flavour g participating with the top in 


16 



























{v,ej) {dpuq) 

Constraint on ebd? 

Observable 

Experimental value 

Viedu 

1.6 X 10 "® 

Rtv 

( 1.230 ± 0 . 004 ) X 10 -'‘ 

Vicdc 

1.4 X 10 "® 

r{D+ -Y ev{) 

< 8.8 X 10 -® 

ViCSU 

1.5 X 10 -® 

BR{K+ -p- evi) 

( 1.55 ± 0 . 07 ) X 10 -® 

UiCSC 

3.9 X 10 "® 

BR{Df — >• ePi) 

< 8.3 X 10 -® 

viebu 

7.8 X 10 -® 

BR{B+ ev) 

< 9.8 X 10 -'^ 

Vicbc 




Vipidu 

3.2 X 10 "® 

Rtv 

( 1.230 ± 0 . 004 ) X 10 -* 

Vipidc 

3.7 X 10 -® 

BR{D'^ — >■ fiv) 

( 3.82 ± 0 . 33 ) X 10-4 

Uifisu 

3.0 X 10 "® 

Rk 

( 2.44 ± 0 . 11 ) X 10 -® 

UifJtSC 

1.5 X 10 -^ 

BR{Df — >■ ftp) 

5.56 X 10 “® 

Vi^bu 

7.7 X 10 -® 

BR{B~^ — ftp) 

< 1.0 X 10 -® 

Pifibc 




ViTdu 

8.0 X 10 "^ 

BR{t — >■ TT+Z/) 

( 1.083 ± 0 . 006 ) X 10 -^ 

ViTdc 

5.2 X 10 "^ 

BR{D+ fp) 

< 1.2 X 10 "® 

ViTSU 

2.6 X 10 -^ 

BR{t K+p) 

( 7.0 ± 0 . 1 ) X 10 -® 

UiTSC 

4.9 X 10 -^ 

BR{D+ fp) 

5.54 X 10 “^ 

ViPbu 

7.2 X 10 -^ 

BR{B^ -f tPt) 

( 1.1 ± 0 . 3 ) X 10-4 

ViTbc 





Table 8: Constraints from “charged current” processes on S' ± P operators. This bound applies to , in an 
basis using tensor operators. The first column is the index combination ijPq^ the second is the constraints, which 
arise from the observable given in column 3. The experimental value used is the last column. Vi is any flavour of 
neutrino. 


the contact interaction. It is given in Table 9. It is clear that the loop gives significant suppression, due to CKM 
elements and light quark masses. 

The VA operators , which involve doublet top quarks, have a larger mixing to down-type quark operators. 
Dressing with a W loop as in Figure 1 gives the operator = (Ji^a Pyij) {daj°‘PLdp), with a log divergent 

coefficient, unsuppressed by light quark masses: 


iia/3 


^YL — 


327r2 


( loc 

r 2 1 

- log 

FA^I 


1 2 2 
\mf - 


m^ 

) 


^YL — ■ 


3Qie 


47r 


-VtpV* 


Jjqt 


(48) 


where (/?, a) € {(6, s), (6, d), (s, d)}, q € {c, u}, only terms with logs were retained, and A ~ Suit in the log to obtain 
the second approximation® and the numerical factors in Table 9. The top operators can also generate Charged 
Current (CC) operators {Vi'^^PL^j) {uq^^PLdp), suppressed only by Vtp, however the bounds on the coefficients of 
such CC operators are weaker than the bounds on FCNC operators, so the best limit on the top operators arises 
from the loop given in Eqs. (47,48). 

The S ± P operators involving either (X = L) or ta {X = R), will generate S ± P operators for 

down-type quarks = (ciPYej) {daPndp), where it is always the heavier dp that is right-handed, because 

there is a mass insertion on an external quark leg. The lepton current is unaffacted by the loop, so will have the 
same chirality in the induced down-quark operator, as it had in the top-LFV operator. There is also a mass insertion 
on the right-handed internal quark line, so the coefficient is oc mt in the presence of ’ : 


ijaP 

^S±P,R 


g'^mtmpVtpV*^ 


qoi 


327r^(m^ — 


1 - 


‘‘W 


" 1 ? - '^W 


log- 


m 


VK 


Stt 


■VtpV^ 


:ic ijqt 


mt 


^S±P,R 


ijqt 

^S±P,R 


LGR LGR 

The result in the presence of is similar: the LQ operator for down-type quarks, 

is generated with a coefficient obtained by replacing eg|.p^ — r ^ 9 - ( 49 )- 


(49) 


{doiPRCj^ 


^We will obtain a bound of e 0.01, which corresponds to A ^ lOmt. We nonetheless conservatively take A 3mt to allow for New 
Particle masses mjvp that are lower than A oc mjv p / Qnp • 
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In the presence of (or the loop factor is oc rugrrip, which is negligibly small. So we obtain no 

interesting bounds on the operators. However, the operators (and contain II and clj, 

between whom the W could be exchanged as indicated in the right diagram of Figure 1. This dressing of the top- 
operator will generate and a tensor operator, in the linear combination corresponding to of 

Equation (21). The loop coefhcient is unsuppressed by masses, as in Equation (48): 


^CC,LQ 






( lo- 

■ m2 ■ 

- log 

lA^I 


I 2 2 

\mt - 



) 


Jjqt 

^S+P,L 


2 a 


em T ^ ijqt 


tl3^S+P,L I 


^CC — 


Vtpe 


ijqt 

LQ,L 


Or, in a basis with tensor operators 


cc 


327r2 


( lo- 

■ ■ 

- log 

FA^I 


1 2 2 



) 


e 


S+P,L — 


TT 


Vtpe 


ijqt 

S+P,L 


(50) 


(51) 


top op. 

down coeff. 

q = u 

q = c 

qXR 

^ijqt 

ij db 
^XL 
ijsb 
^XL 
ijds 
^XL 

8.0 X 10“® 
1.8 X 10“® 
3.2 X 10-® 

7.7 X 10-® 
3.4 X 10-® 
3.1 X 10-^ 

qXL 

^ijqt 

ij db 
^XL 
ijsb 
^XL 
ijds 
^XL 

2.0 X 10“® 
4.6 X 10-4 
8.2 X 10“® 

4.6 X 10-4 
2.0 X 10-® 
1.9 X 10-® 

^^S±P,R 

'^ijqt 

ij db 

^S±P,R 

ijsb 

ijds 

^S±P,R 

2.4 X 10“® 
5.3 X 10-® 
4.1 X 10"® 

5.3 X 10-® 

2.4 X 10-® 
1.8 X 10-® 

.^S+P,L 

'•^ijqt 

®cc 

Mqd 

®cc 

2.5 X 10“® 
1.0 X 10-4 
2.0 X 10-® 

2.5 X 10-® 
1.0 X 10-4 
2.0 X 10-® 


Table 9: The last two columns give the numerical value of the coefficient given in the second column, generated 
by loop containing the top LEV operator given in the first column, of coefficient 2^/20p- The top LEV operator 
involves the top and and the quark Uq G {u, c}; the two choices of q identify the last two columns. The operators 
Ofjqf'^ also induce FCNC operators (involving two down-type quarks), but the coefficients are not listed because 
they are suppressed by two light quark masses. 


C.3 Bounds and cancellations 

The numbers in Table 9 can be compared to the bounds on the various es for light quark operators, given in the 
Tables 6, 7 and 8. For instance, the last column of the third line of Table 9 says that ix%^^(JilaPx^j)(ci'^Ppt) 

generates Pps) with 

= 3.1 X ^ 3.0 X 10-^ ^ eTn ^ 1 

where ^ 3.0 x lO”"^ is the upper bound from K ^ eft given in Table 6, which sets the constraint ^ 1. 
This bound on the efict coefficient of operators appears in the second row and last two columns of Table 1. 
In the case of erct and firct indices, the best bound on the coefficient of the operators arises from B decays. 
This can be seen because Table 1 gives the origin of the bound as bound on these coefficients from 

Table 6 is from B decays®. 

The bounds of Table 1 implicitly assume that only one LFV top operator is present at a time. This assumption 
depends on the choice of operator basis, so it is interesting to consider the possibility of cancellations among the 

®We obtain < 2 X 10® from the limit on but do not include thes numbers in Table 1 because they are to weak to 

be meaningful. 
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V iL A operators. Many of the bounds arise from pseudoscalar down-type meson decays, which are mediated by 
an operator of the form in which case the bound applies to the combinations 

Since the loop that transforms top-operators into down-type quark operators gives different suppression factors for 
operators involving or we neglect possible cancellations among XL and XR operators. So we neglect the 
possbility of cancellations for V ± A operators. 

In obtaining bounds from (pseudoscalar) meson decays, an operator basis that includes tensors is convenient, 
because the tensors do not contribute [49] to the decays M —>■ liij, where M is the pseudoscalar meson made of 
q\q 2 and ti are neutrinos or charged leptons. So we quote low-energy bounds in a basis with tensors^. 

D Appendix: single t production at HERA 

HERA was an e^p collider with ^/s = 319 GeV, where the HI [30] and ZEUS [31] experiments searched for single top 
production: a{e^p —>■ e^tX). We follow an HI analysis [30], which is outlined in Section D.l. Then in Section D.2, 
we discuss how to translate the results of this analysis to a bound on the LEV process a{e^p —>■ pAtX), and in 
Section D.3, we estimate a{e^p —>■ p^tX) in the presence of the operators of Appendix A. 

D.l The HI analysis 

The HI collaboration set a bound [30]: 

2 3 X 10“^ 

(T(e^p —>■ -|-A) < 0.30 pb =-^- at 95%C'L. (52) 

mi 

using 474 pb~^ of data, and requiring the top decay t —> p^vb ®. HI does not give the luminosity in e^p and 
e“p, but from [29,50], we interpret > 228 pb“^ of e+p data and > 205 pb“^ of e~p data at ^/s = 319 GeV, so we 
approximate that half the luminosity was in e“''p, and half in e~p. 

We would like to use this limit to set a bound on a{e^p —>■ p'^t -I- A). This is possible because HI searched for 
leptonic decays of the top, t —>• p'^vb, and did not require to observe the e^. As a first step in a cut-based analysis, 
HI required isolated muons with px > 10 GeV, and missing energy Lfx > 12 GeV; they found 14 events where they 
expected 11 ± 1.8 Standard Model background. Their signal efficiency, for a magnetic moment ufy coupling, was 
44%. A —>■ p^t interaction mediated by our top operators should pass this first set of cuts, so we would like to 

obtain a limit on the coefficients of our operators from this stage of the analysis. Then HI required the muon to 
be positively charged, and in combination with the jet and Lfx , to reconstruct to a top. Only four events remain 
after imposing these cuts, with 2.1 ± 0.3 expected in the Standard Model. Hi’s signal efficiency was 36% and at 
this point they obtain the limit a{e^p —>■ [p'^Lfr 6)e“'"A) < 0.30 pb. 

D.2 Translating to contact interactions 

To translate the HI bound to our contact interactions, we must address various issues. 

1. We estimate that the bound on a{e^p —>■ {p^l^x b)e'^X) after the first cut, is comparable to the final bound 
obtained by HI, so we use the limit of Eq. (52). To obtain this curious approximation, we estimate the t 
production cross section that could be excluded at 95% C.L. using Poisson statistics for 4 observed events 
with 2 background expected and then we estimate the t production cross section that would give a gaussian 
distribution of events centered 4cr above 11, that is at 18. And we find these cross sections are comparable. 

2. Different operators contribute to e“'"p —^ A and e~p —>■ p~t + X. Since the luminosity is approximately 

equally distributed between e+p and e~p, we suppose that the bound for an individual operator is cr ^ 2 x 0.30 
pb. 

^Alternatively, one can impose pseudoscalar decay bounds on ~ From Eq. (49) and after, the top-operators 

qS+p,R qLQ,r induce respectively S + P and LQ operators involving down-type quarks. So in Table 2 are quoted bounds on 
The tx operators and respectively induce the and operators (see Eq. (50)), which 

mediate the leptonic decays of charged mesons, such as > ue'^. So in this case the bound applies to ris quoted 

in Table 2. 

obtained a more restrictive bound, cr < 0.25 pb, by combining the various top decay modes 
®That is, the number of signal events s that give 0.05 = P{s + b)/P(b) where b = background expected. 
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3. In the case of our operators, a central is produced by the contact interaction, so the top can decay to any 
lepton (we do not allow for hadronic t decays, because they are unlikely to produce I^t > 10 GeV). This 
means the upper bound in the cross section is reduced by 1/3, because HI required that the top decay to 

4. We need an estimate for the efficiency for our operators. HI gives a signal efficiency of 44% after their first cut, 
in part because they detect muons as a pr imbalance in the calorimeter, with a 60% efficiency for detecting a 
muon of pt = 10 GeV. Since our event has a central muon, like Hi’s signal, we take the efficiencies to be the 
same, and assume this makes our bounds uncertain by a factor ~ 2 . 

So in summary, we multiply the HI bound by 2/3, and impose 

1 K X 

cr(e"*"p —>■+ V) < 0.20 pb = —- ^ - at 95%C'L. (53) 

rrit 


D.3 (T{e^p —)■ + X) in the presence of top operators 

The operators of Section A induce a differential partonic cross section for e^u —>■ that does not diverge for 

emission at zero angle. We therefore integrate the partonic cross section over all angles (neglecting experimental 
cuts), to obtain: 




I 12(1 

pxx - tt :- 3 - s 




167rm/ 

, (l-m^/s)^ 

487rm/ 


"h \ 

(.+V) 


\£S-P,X 

( 1 -m^/s)^ 

1927rmf 


2 (l-m 2 /s)= 


mr 


1927rm/ ^ 2 ^ 


Oll , Orr 

Olr , Orl 
Os-p,x 


\e-s+p,x + iLQ,x\^][s + + Re{e5+p,xeLQ.x}[s - m^] 


(54) 

(55) 


Os-p,x , Olq ,x 


To obtain the cross section on the proton, we define x = s/s, and obtain, integrating over CTEQIO [51] parton 
distribution functions: 


[ dxx{l-.295/xffuix) =0.0046 

, f dxxfu{x) = 0.041 

/ 0.33 

J 0.33 

/*! 

/■i 

/ dx{l - .295/xffu{x) = 0.0092 

^0.33 

, / dxfu{x) = 0.096 

J 0.33 


(56) 


where the lower x integration limit is estimated to produce, in the centre-of-mass frame, a top of negligible momen¬ 
tum, and a muon of px > 10 GeV. The (1 — .295/a;)^ represents the kinematic suppression factor (1 — m^/sY- 
The resulting cross sections are (for s/m^ = 3.5): 


(j(e^p —>■ p^t + X) 


0.0046s 


|exxr 

16X772-/ 


ep,ut fieut ^ ryr) 

^XX 1 ^xx ^ 


|£xy| 

48X777/ 


0.0046(s -b m/) ^ e/^yXxY < 0-33 


XY ’ 

efiut /leut 


!S^0.0046(, + m>) <0-66 

l«+F,xl“ + ki«x£o.0046(s + m») ^ • 'wA - 


192X777/ 


(57) 

(58) 

(59) 


fieut 

^LQ,X 


< 0.66 


where X,Y G {L, R} and the bounds are estimated by imposing the bound of Eq. (53). 


D.4 What about at HERA? 

The ^reut poorly constrained by rare decays, whereas e^u —>■ could have occured at HERA. It is 

unclear to the authors how HI would have treated a hadronic r decay, so we conservatively restrict to r —>■ evD 
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decays. Allowing any decay of the top, we use the combined HI bound 

1 q X 

a{e'^p ^ e~^t + X) < 0.25 pb = —-^- at 95%C'L, (60) 

mt 

multiplied by 2/BR{t —>■ evV) ~ 12 , where the 2 is from point 2 of Section D.2 above. So a{e^p —>■ + X) ^ 

3.0 pb, and the approximate bounds on are given in Table 1. 

E Appendix: The CMS search for t —> Zq 

The CMS collaboration sets the limit BR{t —>■ Zq) < 6 • 10“^ [11], by searching for ti production, with one top 
decaying leptonically, and the other decaying to Zq followed hy Z ^ ii. In this Section, we emulate in our framework 
this published CMS analysis, which searches for a final state similar to ours. 

We hence implement in our framework the selection described in [11]. The dehnition of good jets and isolated 
leptons is the same as in our analysis described in Section 3.3. We select events containing exactly 3 isolated charged 
leptons (electrons or muons), of which two have the same flavor and opposite sign. Events are requested to contain 
at least 2 jets, and exactly one b-tagged jet, and the missing transverse energy has to be higher than 30 GeV. The 
invariant mass of the opposite sign and same flavor lepton is required to lie within 78 and 102 GeV/c^. 

The two opposite sign and same flavor leptons are identified as coming from the Z decay, if there is more than 
one possible lepton pair, the one with invariant mass closest to the Z mass (mz) is chosen. The remaining lepton 
is associated with the decay of the W and used, with the missing transverse energy, to calculate the neutrino 
longitudinal momentum with a W mass constraint, as explained in Section 3.3. The invariant mass of the W and 
the b-tagged jet is required to be within 35 GeV/c^ of the nominal top mass. The 4-momenta of the two leptons 
coming from the Z are combined with all (non b-tagged) jets in the event. In order for the event to be selected, at 
least one of these combinations must have a resulting invariant mass within 25 GeV/c^ of the nominal top mass. 

All the background processes taken into account in the CMS analysis are considered. In our simulated samples, 
the expected number of background events after applying the CMS selection is 1.5 ± 0.2, to compare with the 
expectation from simulation of the public CMS analysis, of 3.2± I.2(stat)± 1.5(syst). This comparison is not meant 
to be rigorous, as the Delphes emulation of the CMS detector reconstruction is known to be imperfect, and the 
samples used here do not include any effect from pileup. Nevertheless, the agreement is good enough to validate 
our framework as a tool to extract reasonable sensitivity studies. The result of this comparison also motivates the 
variation of 100% in the expected background events for the limit calculation, performed in the main study (Section 
3.4). 

The limit setting procedure is validated as well: when using the number of expected events from the public 
result in our statistical procedure, we obtain an expected upper limit of BR(t —>■ Zg) < 0.9 • 10“^ to compare to 
the public result BR{t —>■ Zq) < 1 • 10“^. 

Having validated our method on expected background allows us to evaluate the constraint set on the t qp^e^ 
branching ratio by the CMS analysis [11]. Indeed, our LEV top events and the CMS analysis share the same overall 
final state (missing transverse energy, 3 isolated leptons, and 2 jets, of which one is a 6-jet). Half of our LEV events 
should give a pair of opposite sign, same flavour leptons, as required in the CMS analysis. However, in our case, 
they do not come from the Z, but respectively from the t and i (see Eigure 2), so they are easily rejected by the 
CMS requirement that their invariant mass be near mz- In order to quantify this, we evaluate the efficiency of the 
CMS event selection on the t —>■ qp^e^ signal sample, obtaining e = 0.050 ± 0.005%. This would correspond to 
an expected limit on LEV in top decays of BR{t —^ qp^e^) < 3.7 • 10“^, according to the limit setting procedure 
based on RooStats (a simple rescaling of the published limit with selection efficiencies and cross sections would give 
a comparable value of BR{t —>■ qp^e^) < 4 • 10“^). Comparing to Equation (4), we see that this is on the verge of 
probing LEV top decays mediated by a four-fermion operator. 
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